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The quantum statistical interference between identical particles emitted from a com- 
pletely chaotic source is expected to provide valuable input for the space time description 
of the system. Intensity interferometry of thermal photons produced in heavy ion col- 
lisions is a very promising tool to explore the structure and dynamics of the collision 
fireball. Thermal photons having Kt < 2 GeV/c get competing contribution from both 
hadronic and quark matter phases and this competition gives rise to a rich structure in 
the outward correlation function, owing to the interference between the photons from the 
two sources. The temporal separation between the two sources provides the lifetime of 
the system and the correlation results are found to be sensitive to quark hadron phase 
transition temperature and the formation time of the plasma. The outward correlation 
function strongly depends on the equation of state of the strongly interacting matter and 
is seen to clearly distinguish between the lattice based and bag model equations of state. 
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I. INTRODUCTION 

The primary motivation of the relativistic 
heavy ion collision studies is to investigate 
the properties of a deconfined strongly inter- 
acting matter, Quark Gluon Plasma (QGP) 
produced in the collisions. The observation 
of elliptic flow [J H|, jet quenching [!, B, 
and the success of parton recombination Q 
as a model for hadronization at the Relativis- 
tic Heavy Ion Collider (RHIC) have advanced 
our understanding of this hot and dense novel 
state of matter (Fig. [T]). The evolution of the 
system produced in the collision is explained 
quite successfully using the powerful method 
of hydrodynamics |6|. However, numerous im- 
portant questions like; how quickly (if at all) 
does the plasma thermalize, what is the life 
time of the system, if there is a phase transi- 
tion or not, are still to be addressed. 

The small size (~fm) and transient (~fm/c) 
nature of the collision fireball make it ex- 
tremely difficult to obtain space time infor- 
mation of the system. Correlation between 
two final state particles which is stronger for 
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smaller separation in space time as well at 
smaller relative momentum, provides the most 
direct tool to get spatio-temporal information 
of the system 

In 1950s', Robert Hanbury Brown and 
Richard Q. Twiss measured the angular size of 
a distant star by utilizing the method of inten- 
sity interferometry for the very first time Q. 
Theoretically the photon bunching was first 
explained by Purccll [9] , in one of the key ex- 
periments of quantum optics. In 1959 Gold- 
haber et al. [lOj observed an unexpected an- 
gular correlation among identical pions while 
studying the po resonance, which was the first 
observation of intensity correlation in particle 
physics. Very soon it was realized that the 
correlations of identical particles emitted by 
highly excited sources are sensitive not only 
to the size of the system, but also to its life- 
time and the momentum dependent correla- 
tion function contains information about the 
dynamics [ll[ of the system produced in the 
collisions. 

The heavy ion community often uses the 
term HBT in reference to the original work 
of Hanbury Brown and Twiss for any analy- 
sis related to the spatio-temporal aspects of 
the particle emitting source. Several interest- 
ing theoretical works on correlation study in 
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FIG. 1: Elliptic flow (upper panel) PJ], jet 
quenching (middle panel) [l5j, and the success of 
parton recombination as a model for hadroniza- 
tion (lower panel) [l(J at the RHIC have advanced 
our understanding of the hot and dense medium 
produced in relativistic heavy ion collisions. 



heavy ion collisions were reported in 1970s' by 
Shuryak [l2j], Gyulassy et at [l3[ and many 
others. Other important contributions in the 
eighties include parametrization of the source 
function and a more detailed analysis of the 
role of final state interactions in the system 
(see review article Q for detail). 
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FIG. 2: Thermal photon predictions for central 
Au+Au reactions at RHIC computed with differ- 
ent hydrodynamical and dynamical fireball mod- 
els along PHENIX experimental data fl7| ]. 



II. INTENSITY INTERFEROMETRY 
OF PHOTONS 



Electromagnetic radiation, known as the 
thermometer of the medium, has long been 
considered as one of the most promising and 
efficient probes for the hot and dense state 
of matter produced in the relativistic collision 
of heavy nuclei []~7j j . Intensity interferometry 
with identical hadrons or light nuclei has been 
an important tool to learn about the dynamics 
of the subatomic or nuclear collisions. How- 
ever, the hadrons suffer final state interactions 
and their correlations mainly carry informa- 
tion about the late dilute stage of the system. 

Photons on the other hand, decouple from 
the system immediately after their creation, 
suffering negligible re-scatterings with the 
medium (a <C a s ), and carry undistorted 
information about the circumstances of their 
production to the detector. Obtaining direct 
information about the earliest hot and dense 
stage of the system by studying large Kt pho- 
tons is very promising. A large production 
of photons having very high transverse mo- 
menta is expected from the pre-equilibrium 
stage, and some additional sources of photons 
have also been proposed recently [H, [h| . The 
major problem of-course arises from the mea- 
ger emission of direct photons compared to 



the huge background of decay photons (mostly 
from 7T° and ij mesons) produced in the colli- 
sions. Of late, there have also been tremen- 
dous advances in methods for identification of 
single photons [20]. 

The theory of the intensity interferometry 
of photons from relativistic heavy ion colli- 
sions has been pursued in considerable detail 
by several authors [2ll - l23| in last two decades. 
These studies further refined the early ex- 
pectations of photon as powerful tool to get 
spatio-temporal aspects of the fireball evolu- 
tion. In a recent study by Frodermann and 
Heinz (24[, the details of angular dependent 
intensity correlation of thermal photons from 
non-central collision of heavy nuclei as well 
the dual nature of thermal photon emission 
resulting from the superposition of QGP and 
hadron resonance gas photon production have 
been highlighted. 

It is generally felt that the experimental 
efforts for these studies have a larger likeli- 
hood of success at RHIC and LHC energies 
because of larger initial temperature of the 
plasma and a large suppression of pions due 
to jet- quenching. 

Till date only one measurement of photon 
intensity correlation at very low transverse 
momentum Kt has been reported by WA98 
collaboration for central collision of lead nu- 
clei at CERN SPS 251. 



of photons, EdN/d 4 xd 3 k from the quark and 
hadronic matter phases produced in the col- 
lisions and the correlation function C(q, K) 
is decomposed in terms of the outward, side- 
ward, and longitudinal momentum differences; 
q , q s and qi respectively. One can write the 
four-momentum of the ith photon in-terms of 
transverse momentum kx, rapidity y, and az- 
imuthal angle ip as, 

k ; = (k iT cos ipi, k iT sin fy, k iT sinh yi), (2) 

and the three momentum differences q a , q s and 
qe are of the form: 
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— k\T sinh yi — k 2 T sinh y 2 ■ (3) 

The radii corresponding to these momentum 
differences are obtained by approximating the 
correlation function to Gaussian parametriza- 
tion as; 



A. Formulation 

The spin averaged intensity correlation be- 
tween two photons with momenta ki and k2, 
emitted from a completely chaotic source can 
be expressed in terms of the space-time emis- 
sion function S(x,K.) as (2ll. [2q: 



1 



C(q,K) = l+- 



|J d 4 x S(x, K)e i 



2jd i xS(x,k 1 ) Jd*xS(x,-k 2 ) 

(1) 

where the obtained correlation is a function of 
the relative q(= ki — k2) and average mo- 
menta K (= (ki + k2)/2) of the two pho- 
tons. The space-time emission function S is 
often approximated as the rate of production 



C(q , q s , qe) = 1 + \ cxp [- (q 2 R 2 + q 2 R 2 s 

+ 8?ii?)].(4) 

The root mean square momentum difference 
(qf ) and the radii are obtained as, 



J(C-l)qfd qi 
J(C-l)dq % 
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It is to be noted that the l/[2 (<jf )] 1 ' 2 becomes 
a useful measure when the correlation function 
has a more complex nature. 
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FIG. 3: Outward (top panel), sideward (middle 
panel), and the longitudinal (lower panel) corre- 
lation functions for thermal photons produced in 
central collision of gold nuclei at RHIC. Symbols 
denote results of the correlation; curves denote 
fits [H. 



B. Initialization and system evolution 

The exciting possibility of observation of the 
interference of thermal photons in the inter- 
mediate Kt (< 2 GeV/c) range has been ex- 
plored in detail in a recent work of Srivastava 
et al. [26]. Photons in this momentum range 
are unique as they get competing contribution 
from the two phases; they have their origin ei- 
ther in the hot and dense quark phase of the 
system or in the relatively cooler but rapidly 
expanding hadronic phase, where a large build 
up of the radial flow boosts their transverse 
momentum. This competition gives rise to a 
rich structure specially in the outward corre- 
lation function, owing to the interference be- 
tween the photons from the two sources. 

For this particular study, central collisions 
of gold and lead nuclei at the top RHIC and 
LHC energies respectively are considered. It 



is assumed that a thermally and chemically 
equilibrated plasma is formed at a very small 
initial time tq (0.2 fm/c for RHIC and 0.1 
fm/c for LHC). An isentropic expansion of the 
plasma is considered where the initial temper- 
ature Tq is obtained from the relations, 



2tt 4 1 dN 



4aT 3 r 



(7) 



here At is the transverse area of the sys- 
tem, dN/dy is the particle rapidity density, 
and a = 42.257r 2 /90 for a plasma of mass- 
less quarks and gluons. The initial energy 
density is estimated by combining 25% con- 
tribution from binary collisions and 75% from 
number of wounded nucleons. The bag model 
equation of state is considered for the sys- 
tem evolution where a first order phase tran- 
sition takes place at a temperature of about 
180 MeV and the freeeze-out temperature is 
taken as 100 MeV. The relevant hydrodynamic 
equations are solved under the assumption of 
boost invariant longitudinal and azimuthally 
symmetric transverse expansion using the pro- 
cedure discussed in Ref. 27J and integration 
performed over the history of evolution. Fig. [2] 
shows thermal photon results from from hy- 
drodynamic model by different groups along 
with direct photon data from PHENIX. 

The complete leading order results for pro- 
duction of photons from quark matter by 
Arnold et al. [28| , and the results of Turbide et 
al. [29I] for radiation from the hadronic matter 
are used for this analysis. The values of charge 
particle multiplicity at mid rapidity are taken 
as 1260 [HI for 200 A GeV Au+Au collisions 
at RHIC and 5625 [H for 5.5A TeV Pb+Pb 
collisions at LHC. 



C. Correlation functions for thermal 
photons 

The results for the outward, sideward, and 
longitudinal correlation functions for thermal 
photons at RHIC having Kt ~ 2 GeV/c are 
shown in Fig. |3] The four-momenta of the two 
photons are chosen so that when the outward 
correlation is studied, q s and qe arc identically 




FIG. 4: Transverse momentum dependence of the 
outward radii and temporal duration obtained by 
fitting the final correlation function for thermal 
photons at RHIC. Radii determined from he root- 
mean-square momentum difference for the corre- 
lation function are also given for comparison [26|. 

zero and the dependence on qo is clearly seen, 
and so on. Considering only the quark matter 
or the hadronic matter only, it is found that 
the correlation functions for the two phases 
can be approximated as 

C{ qi ,a) = 1 + 0.5\ Pha \ 2 (8) 

where i = o, s, or £, and a denotes quark 
matter (Q) and hadronic matter (H) in an ob- 
vious notation. The source distributions \pt, a \ 
is well described as a Gaussian function, 

|pi, a | = I a exp [- 0.5 (g?i$ J] (9) 

where I Q = dN Q /(dN Q + dN H ) and I H = 
dNn/{dNQ + dNjj) are the fractions of the 
photons from quark matter and hadronic mat- 
ter respectively. 

The final correlation functions are then ap- 
proximated as, 

C{ qi ) = 1 + 0.5 [K Q | 2 + K ff | 2 
+ 2\p itQ \\p itH \ cos(qiARi)} 

(10) 

where the cosine term in the above equation 
clearly brin gs o ut the interference between the 
two sources [31| ■ Here ARi stands for the sep- 
aration of the two sources in space and time 
and q is the four momentum difference. For 



FIG. 5: Transverse momentum dependence of the 
fraction of thermal photons from the quark matter 
and hadronic matter phases at RHIC with varying 
transition temperature [26| . 

K T ~ 2 GeV/c, at RHIC, the value of AR 
is 12.3 fm obtained from fitting whereas AR S 
and ARi are of the order of zero. 

These results imply that, while the spa- 
tial separation of the two sources is negligi- 
ble, their temporal separation is much larger, 
which gives the lifetime of the system. If the 
mixed phase is of shorter duration or absent, 
this will obviously decrease. 

The Kt dependence of AR a , and the out- 
ward radii for the hadronic and quark mat- 
ter sources of photons obtained from fitting 
are shown in Fig. [4] which clearly explains the 
fact mentioned above. The outward radius 
for the quark contribution depends weekly on 
the transverse momentum, which is indicative 
of a mild development of the radial flow dur- 
ing the quark matter phase. The correspond- 
ing radius for the hadronic contribution shows 
a stronger dependence on the transverse mo- 
mentum resulting from a more robust devel- 
opment of radial flow during the late hadronic 
phase. 

A very interesting structure is observed for 
the duration of the source, which increases 
slightly at higher Kt as the momenta of the 
photons emitted from the hadronic phase is 
blue shifted to much larger value due to strong 
radial flow. The saturation of AR Q towards 
low Kt has its origin in the competition 
between radial expansion and decoupling of 
hadronic matter as it cools down below the 
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FIG. 7: Tc dependence of the outward correlation 
function at RHIC H. 



temperature of freeze-out at the edges. 

As seen from Eq. ([6]) , the inverse root mean 
square momentum, which is a measure of the 
correlation radius, is seen to vary rapidly with 
transverse momentum for the outward correla- 
tion owing to rapid variation of the competing 
contribution of photons from the two phases. 

As the nature of the overall correlation func- 
tion strongly depends on the relative contri- 
bution from the two phases, the fraction of 
the quark (Iq) and hadronic (Ih) contribu- 
tions are investigated as a function of trans- 
verse momentum as well as varying critical 
temperature (see Fig. [5J . One can see that the 
transverse momentum at which the quark and 
hadronic contributions become equal is sensi- 
tive to the transition temperature and a de- 
creasing Tc increases the fraction of photons 
coming from the quark matter. 

The variation of the source function as a 
function of transverse distance and time for 
different values of transverse momentum has 
been studied in detail (not shown here) in or- 
der to understand the correlations results and 
the relative contributions of the two phases in 
it. It is observed that the radial distribution 
for both the phases are centered at tt near 
zero and the source function for the hadronic 
phase extends considerably beyond the same 
for the QGP phase. This is of course due 
the large transverse expansion of the system 
and results in a good amount of thermal pho- 
ton production from hadronic phase at larger 
radii. A more valuable insight is provided by 



the temporal distribution of the source func- 
tion, where a large duration of the hadronic 
phase is observed which controls the parame- 
ter AR Q . 



D. Sensitivity to formation time and 
transition temperature 

We know that a smaller value of the initial 
formation time of the plasma implies a larger 
initial temperature. The sensitivity of the cor- 
relation functions to the initial formation time 
of the plasma is investigated by considering 
systems with identical entropies but varying 
formation time. Smaller ro leads to a larger 
production of photons from the quark matter 
phase and as a result the outward correlation 
function changes its nature as shown in Fig. |6] 

It is to be noted that a larger value of for- 
mation time (~ 1 fm/c) would necessitate 
inclusion of the pre-equilibrium contribution 
of the photons, that must surely be there at 
least at larger values of the transverse momen- 
tum. Another important contribution, pho- 
tons from jet conversion processes also likely 
to contribute at higher Kt- These contribu- 
tions would increase the fraction of quark mat- 
ter contribution in the correlation measure- 
ment, which in turn could mimic an effective 
smaller value of the formation time. 

In another potentially powerful observation 
the sensitivity of the outward correlation func- 




FIG. 8: Speed of sound for (i) hadron resonance 
gas, (ii) volume corrected hadron resonance gas, 
(iii)hadron and Hagedorn resonance gas, and (iv) 
volume corrected hadron and Hagedorn resonance 
gas with lattice results [32]. 



FIG. 9: Pion pr spectra for the equation of state, 
HHB and HHL for central collisions of Au nuclei 
at RHIC. Experimental data for pions (0-5% cen- 
trality bin) [36|] are shown for comparison [3^ ]. 



tion to the transition temperature Tc is shown 
in Fig. where the change in the interference 
pattern at larger qo can be easily understood 
in terms of the relative life time of the two 
phases. 

The correlation results for Pb+Pb collisions 
at LHC energy also show similar qualitative 
nature as observed at RHIC. As the initial 
temperature likely to be attained at LHC is 
much larger than that at RHIC, the study of 
the properties of QGP phase and the dynam- 
ics of its evolution will be more interesting at 
LHC. A higher initial temperature would lead 
to a longer duration of the interacting system, 
which in turn would provide ample opportu- 
nity for the mechanism of expansion to de- 
velop. 

III. EQUATION OF STATE AND 
INTENSITY CORRELATION OF 
THERMAL PHOTONS 

In a very recent work by De et al. [32| the 
equation of state (EOS) of strongly interact- 
ing matter and how the change in EOS af- 
fects the results of system evolution, particle 
spectra as well as the intensity interferometry 
of thermal photons have been explored in de- 
tail. It is shown that an equation of state for 
hot hadronic matter consisting of all baryons 



having M < 2 GeV and all mesons having 
M < 1.5 GeV, along with Hagedorn reso- 
nances [33l ] in thermal and chemical equilib- 
rium, matches rather smoothly with lattice 
equation of state (p4 action, N T — 8) [34[ 
for zero baryon chemical potential for T up 
to ~ 200 MeV, when corrected for the finite 
volume of hadrons. 

Furthermore, two equations of state for 
strongly interacting matter are constructed; 
one, HHL, in which the above is matched to 
the lattice equation of state at T = 165 MeV 
and the other, HHB, where it is matched to a 
bag model equation of state with critical tem- 
perature T c — 165 MeV. The lattice equa- 
tion of state displays a sharp cross-over for 
180 < T < 190 MeV. 

The results for the energy density, pressure, 
and the entropy density of the hadronic matter 
with successively increasing richness of the de- 
scription; viz., hadron gas, hadron+Hagedorn 
gas, and volume corrected hadron+ Hagedorn 
gas have been studied in detail (see Ref. [32[ 
for detail). Best agreement with the lattice 
calculation is obtained when the hadron + 
Hagedorn gas is corrected for the finite volume 
of the particles and as a result the hadronic 
matter for both HHL and HHB equations of 
state is constructed by combining hadron and 
Hagedorn gas with volume correction. The re- 
sults for the square of speed of sound for the 
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FIG. 10: The outward correlation function for 
thermal photons at RHIC [321 ]. 



four description of the hadronic matter and 
their comparison with the one obtained from 
the lattice calculations are shown in Fig. [8] 



A. Particle spectra for HHB and HHL 
equations of state 

Considering an ideal hydrodynamic (az- 
imuthally symmetric and longitudinally boost 
invariant) evolution of the system and the 
same initial conditions as used in Ref. [26[ 
for Au+Au collisions at RHIC as well as for 
Pb+Pb collisions at LHC, transverse momen- 
tum spectra of various hadrons and photons 
are calculated for the two equations of state. 
The freeze-out is assumed to take place at 
a temperature of about 100 MeV and the 
particle spectra are obtained using Cooper- 
Frye formulation [35[. Whereas for the pho- 
ton production, the standard rates from the 
Refs. jH H<| are used. 

Transverse momentum spectra of pions and 
thermal photons for central collision of Au nu- 
clei at RHIC and at midrapidity are shown in 
Figs. [9] and [TT] along with experimental data 
from PHENIX [H, |37| for the two cases. It 
is observed that both the equations of state 
give a reasonable description of the particle 
distribution with a slight preference for the 
HHL EOS and the inverse slope of the spectra 
for HHL EOS is larger than the same for the 
HHB EOS. Similar results have been reported 




FIG. 11: Thermal photon spectra for the equation 
of state, HHB and HHL for central collisions of 
Au nuclei at RHIC. Experimental data for (0-20% 
centrality bin) [33] are shown for comparison [32] ■ 

for the bag model and lattice EOS in Ref. [38| 
at RHIC earlier. It is to be noted that a com- 
plete description of photon data would involve 
addition of prompt contribution scaled by ap- 
propriate nuclear overlapping function. 

Corresponding results for the LHC (see 
Ref. [32] for detail) show that the difference 
between the transverse momentum distribu- 
tion of photons and the hadrons for the two 
EOS is further reduced. However, a slight in- 
crease in the inverse slope for the HHL EOS 
is observed as a result of the larger lifetime of 
the system at the LHC energy. 

One can conclude from these results that 
particle spectra (even for thermal photons) 
can not distinguish between the equations of 
state, one admitting a first order phase transi- 
tion and the other admitting a rapid cross-over 
as suggested by lattice calculations. However, 
it is to be noted that the history of evolution 
of the two systems is different as they are sub- 
jected to different rates of expansion due to 
the varying speed of sound. 

In the same study it is observed that the re- 
sults for the temporal evolution of average en- 
ergy density, temperature and radial flow ve- 
locity for both the RHIC and LHC energies do 
not show significant difference for the two EOS 
HHL and HHB. The most profound variation 
is observed for the radial velocity for the two 
EOS, where the radial velocity for HHL rises 
continuously, stays below that for the HHB, 
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FIG. 12: Transverse momentum dependence of 
the outward radii for hadronic and quark matter 
sources for thermal photons (32T |. 

but overshoots it once the later stalls due to 
the onset of the mixed phase. As a result, fi- 
nal Vt for the HHL EOS is slightly larger than 
that for the HHB. 



A very clear difference in the underlying 
radii for the hadronic matter contribution 
(shown in Fig. H2|) for the two equations of 
state is observed; where the one for the HHB 
equation of state being effectively much larger. 

One can see that the difference AR Q , which 
corresponds to the lifetime of the system is 
slightly smaller for the HHL equation of state, 
as it does not incorporate mixed phase like the 
HHB. The R for the hadronic matter for the 
HHL equation of state is also much smaller for 
the same reason. 

Similar qualitative nature for the correla- 
tion functions and radii is observed at LHC 
energy also. In view of these important obser- 
vations, it would be quite interesting to ana- 
lyze this behavior using a fully 3+1D hydro- 
dynamics |39j . 

IV. SUMMARY AND CONCLUSIONS 



B. EOS and intensity interferometry of 
thermal photons 

It is observed that the HHL equation state 
leads to a larger production of photons at 
smaller radial distances as well as at inter- 
mediate times. We know that the photons 
from quark matter originate at early times and 
those from hadronic matter are emitted from 
the later part of the system evolution As the 
interference between the two depends on their 
relative contributions, this difference holds out 
a promise that we may see a difference in the 
intensity interferometry of photons for the two 
equations of state. 

It is found that the sideward correlation 
is essentially identical for the two equations 
of state at both RHIC and LHC energies, 
whereas the longitudinal correlation function 
shows a slight variation when the EOS is 
changed from HHB to HHL (see Ref. [H for 
details). The results for the outward corre- 
lation (shown in Fig. ITU1) are more dramatic 
and show a clear difference for the two equa- 
tions of state. The same treatment of Ref. [26| 
is followed to estimate the source size of the 
hadronic and the quark matter contributions. 



The interference of thermal photons from 
the quark matter and hadronic matter phases 
produced in relativistic heavy ion collisions 
gives rise to a unique structure, specially in 
the outward intensity correlation function for 
intermediate transverse momenta. The corre- 
lation results are found to be quite sensitive 
to the initial conditions and the equation of 
state of the strongly interacting matter. This 
can be used to distinguish between the bag 
model and lattice based equations of state. 
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